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A mode-coupling theory for the slow single-particle dynamics in fluids adsorbed in disordered 
porous media is derived, which complements previous work on the collective dynamics [V. Krakovi- 
ack, Phys. Rev. E 75, 03f503 (2007)]. Its equations, like the previous ones, reflect the interplay 
between conflnement-induced relaxation phenomena and glassy dynamics through the presence of 
two contributions in the slow part of the memory kernel, which are linear and quadratic in the den- 
sity correlation functions, respectively. From numerical solutions for two simple models with pure 
hard core interactions, it is shown that two different scenarios result for the diffusion-localization 
transition, depending on the strength of the confinement. For weak confinement, this transition 
is discontinuous and coincides with the ideal glass transition, like in one-component bulk systems, 
while, for strong confinement, it is continuous and occurs before the collective dynamics gets non- 
ergodic. In the latter case, the glass transition manifests itself as a secondary transition, which can 
be either continuous or discontinuous, in the already arrested single-particle dynamics. The main 
features of the anomalous dynamics found in the vicinity of all these transitions are reviewed and 
illustrated with detailed computations. 



I. INTRODUCTION 

The dynamics of fluids confined at the nanoscale is a 
topic of great interest, both for fundamental and prac- 
tical reasons [1, 2, 3]. Among the relevant questions 
from a fundamental point of view, one counts anomalous 
molecular transport characterized by subdifFusive laws, 
which result from the trapping and obstruction phenom- 
ena due to geometric and topological constraints and/or 
quenched disorder [4, 5, 6]. Another topic which has 
attracted considerable attention in the past few years is 
the dynamics of confined glass- forming liquids [7, 8] , with 
the main motivation of clarifying and substantiating the 
concept of cooperativity, a key ingredient of many glass 
transition theories [9, 10, 11]. On the practical side, the 
study of many problems in the fields of biology, geology, 
and chemical engineering could benefit from a finer un- 
derstanding of the dynamical behavior of confined fiuids. 

Recently, as a first step towards a general and versatile 
theoretical approach to the problem, we have developed 
an extension of the ideal mode-coupling theory (MCT) 
of the liquid-glass transition [12, 13, 14, 15, 16] which 
allows one to study the slow collective dynamics of fluids 
adsorbed in disordered porous solids on the basis of the 
model of the so-called " quenched-annealed" (QA) binary 
mixture [17, 18, 19]. In these systems, first introduced 
by Madden and Glandt [20, 21], the porous solid (the 
quenched component) is represented by a random array 
of particles frozen in a disordered configuration sampled 
from a given probability distribution, in which the fiuid 
molecules (the annealed component) equilibrate. 

The theory makes a number of interesting predictions 
[17, 18, 19]. Indeed, two distinct types of ideal liquid- 
glass transition scenarios are predicted, which are either 
discontinuous, like in the bulk, in situations of weak con- 
finement (dilute matrices), or continuous in situations of 
strong confinement (dense matrices). In the intermediate 



region where the nature of the transition changes, degen- 
erate or genuine higher-order singularities and glass-glass 
transition lines are found depending on the details of the 
fluid-solid system. Moreover, a reentrant glass transi- 
tion line is predicted in the low fluid-high matrix density 
regime. 

In the present paper, we address within the same the- 
oretical framework the dynamics of a tagged particle 
moving in a fluid adsorbed in a disordered porous solid. 
More precisely, we shall derive and solve the MCT equa- 
tions describing the time evolution of the wave-vector- 
dependent tagged-particle density correlators and, from 
their small wave vector limit, of the mean-squared dis- 
placement. These quantities can be measured by a vari- 
ety of techniques and usually more easily than the corre- 
sponding collective properties. For instance, because hy- 
drogen has a large incoherent scattering length, inelastic 
neutron scattering techniques essentially probe the spec- 
trum of the incoherent intermediate scattering function 
when applied to hydrogen-rich molecular systems like wa- 
ter and many organic glassformers; deuterated samples 
are required to measure its coherent counterpart. Also, 
in computer simulations, tagged-particle quantities can 
be averaged over all fluid particles in the sample, leading 
to a very significant improvement in statistics compared 
to the collective functions. Thus, we anticipate that with 
the present developments it will become easier to com- 
pare the predictions of the theory with high quality ex- 
perimental and simulation data. 

A priori, one expects to see manifestations of two types 
of ergodicity breaking events in the present problem. The 
ideal liquid-glass transition as obtained from the study of 
the collective dynamics is the first one, following the same 
mechanism as in the bulk: the tagged-particle dynamics 
is coupled to the collective dynamics and a spontaneous 
arrest of the latter might trigger a similar arrest of the 
former. The second one is the diffusion-localization tran- 



sition corresponding to the permanent trapping of the 
tagged particle in some domain of finite spatial extent. 
This is a confinement- and disorder-induced phenomenon 
which is already present when the tagged particle is alone 
in the porous medium. It is one aim of the present work 
to understand how these two processes interplay. Note 
that early MCT studies of bulk binary hard sphere mix- 
tures have clearly shown that the relation between the 
two phenomena can be nontrivial [22, 23, 24]. 

From the point of view of the diffusion-localization 
transition, the reported theory represents an exten- 
sion of early work on the classical random Lorentz gas 
[25, 26, 27], incorporating collective effects due to finite 
fluid densities at two levels. First, through the use of the 
proper statistical mechanical formalism for the definition 
and the computation of the static correlation functions of 
the QA mixture which are needed as input of the dynam- 
ical theory, the dependence of the structure of the fluid 
on its density and on the fluid-fluid interactions is nat- 
urally accounted for. Second, as already mentioned, the 
coupling of the tagged-particle dynamics to the collective 
dynamics is explicitly included. So, in addition to the lo- 
calization effect of the porous medium, obstruction phe- 
nomena due to molecular crowding or, conversely, decor- 
relation processes due to fluid-fluid collisions can in prin- 
ciple be captured by the theory. 

The paper is organized as follows. In Sec. II, the model 
of the QA mixture is defined, the MCT equations for the 
collective dynamics are recalled and those for the tagged- 
particle dynamics are derived and shortly discussed. In 
Sec. Ill, the relation of the theory to previous work and 
a variant are explored, in an attempt to provide guide- 
lines for future studies with other formalisms. The dy- 
namical phase diagrams are considered in more details 
in Sec. IV. Section V is devoted to the presentation of 
the complete relaxation scenarios and Section VI to con- 
cluding remarks. A few technical results are presented in 
appendices. 



II. THE MODEL AND ITS MODE-COUPLING 
THEORY 

As in MCT studies of bulk systems, the variables of 
interest in the present theory are the Fourier components 
of the microscopic densities characterizing the system. 
For the problem of the dynamics of a tagged particle 
moving in a fluid adsorbed in a disordered porous solid 
(see Fig. 1 for a sketch of the system), there are four such 
quantities which are relevant either to the statics or the 
dynamics. 

In a QA system, the disordered porous medium is rep- 
resented by a collection of A^m rigorously immobile point 
particles, randomly placed in a volume V at positions de- 
noted by si, S2, . . . , sjv„ according to a given probability 
distribution ^(si, S2, . . . , sjv^^J [20, 21]. Its overall den- 
sity is Tim = Nm/V and its frozen density fluctuations 
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FIG. 1: Sketch of a QA system. In black, the immobile ma- 
trix particles. In white and gray, with arrows symbolizing 
their movement, the fluid particles and the tagged particle, 
respectively. 

are given by 

^m^ggiqs,^ (1) 

where q denotes the wave vector. 

The fluid component consists of Ni point particles 
(density rif = Nf/V) of mass m, which equilibrate at 
a temperature T in the random potential energy land- 
scape created by the frozen matrix particles. As usual, 
its time-dependent density fluctuations are defined as 

p^(i) = Ee^'>-^(*), (2) 

where rj{t) is the position of the fluid particle j, j = 
l,...,Nf, at time t. But, in the case of a QA sys- 
tem, these qiiantities arc not immediately useful. In- 
deed, because of the random external field provided by 
the fixed matrix, nonzero average density fiuctuations ex- 
ist at equilibrium. Thus, one is led to consider relaxing 
and non-relaxing fluid density fluctuations, correspond- 
ing to i5/0q(t) = Pq{t) — {Pq) and (/3q), respcctively, where 
(• • • ) denotes a thermal average taken for a given realiza- 
tion of the matrix. 

A tagged particle of mass mg is immersed in the sys- 
tem. Its time-dependent density fluctuations are given 

by 

p^(t)=e^ir=«, (3) 

where r^{t) is its position at time t. In an inflnite system, 
(Pq) ~ so there is no need to separate relaxing and 
non-relaxing parts in this quantity. 

The theory deals with disorder-averaged quantities un- 
der the assumption that the matrix is statistically ho- 
mogeneous, so that, while for any single realization, the 
system lacks translational and rotational invariance, all 
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expectation values computed with the matrix probabil- 
ity distribution have the same properties as in a truly 
translationally and rotationally invariant system. For in- 
stance, (pfj) = 0, where denotes the average over the 
matrix realizations performed after the thermal average 
(•••). It also results that the disorder-averaged correla- 
tion functions of all the density fluctuations present in 
the problem are diagonal in q and only depend on its 
modulus q. 

The collective dynamics has been considered in 
Refs. [17, 18, 19]. It is described by a closed set of self- 
consistent equations for the time evolution of the nor- 
malized connected autocorrelation function of the fluid 
density fluctuations 
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where Sq^'^^ is the connected fluid-fluid structure factor 
deflned as 



(5) 



They consist of a standard generalized Langevin equa- 
tion. 
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dTmq{t - T)^q{T) = 0, (6) 



with il^ = q^ksT /mS^'"'^^ and initial conditions 0<j(O) = 

1, 0g(O) = 0, and of a mode-coupling approximation for 

the memory kernel, mq{t) = Tq5{t) + m^'^\t), where 
Tq is a friction coeflacient associated with fast dynamical 
processes and 
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(1). 
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static correlation functions used in this work are reported 
in Appendix A. 

Wo now turn to the tagged-particle dynamics as en- 
coded by its density correlator 
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Using projection operator methods, a generalized 
Langevin equation is first obtained, 
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(/)^(0) = 0. This equation is exactly the same as for bulk 
systems. 

It now remains to derive a mode-coupling approxima- 
tion for the memory kernel m^(f). Here, some care is 
needed. Indeed, it was found in Ref. [19] that, in the 
case of an isolated molecule moving in a random porous 
matrix, a direct calculation based on the fact that the 
molecule interacts with the solid only [25, 26, 27] does 
not yield the same result as taking the limit of vanishing 
fluid density in the equations describing the collective dy- 
namics, while the two approaches should be equivalent. 
The present problem, of which the above is actually a spe- 
cial case, shows a similar difficulty. A direct approach in 
which the tagged-particle density fluctuations only cou- 
ple to the matrix and fluid density fluctuations results 
in different equations as considering an adsorbed binary 
fluid mixture in the limit of a vanishing concentration of 
one of the components. In the absence of any obvious 
physical reason to consider the zero concentration case 
as a special one, we make the choice of consistency at all 
densities and use the latter strategy. Additional argu- 
ments in favor of this choice and the equations obtained 
with the alternative approach are reported in the next 
section. 

So we start with the MCT equations for a mixture 
adsorbed in a disordered porous matrix, which are given 

)r completeness in appendix B, specialize them to the 
,se of a binary mixture, and send the concentration of 
le component to zero. Keeping only the leading terms, 
we get mq(f) = T\6{t) + mq^^'^\t), where here again 
is a friction coefficient associated with fast dynamical 
processes and 
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Sq^'"^ is the disconnected or blocked fluid-fluid structure 
factor 
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and Cg is the Fourier transform of the connected fluid- 
fluid direct correlation function [28, 29, 30, 31]. Its 

blocked counterpart will be denoted c, 
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in the follow- 



ing. For reference, the relations between the different 
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Cq''' and hq^'^^^ arc the Fourier transforms of the con- 
nected singlc-particle-fluid direct correlation function 
and of the single-particle-single-particle blocked total 
correlation function, respectively (see Ref. [32] for general 
definitions). They obviously have blocked and connected 

counterparts, denoted by c^q^^^ and h^q^''\ respectively. 

Like in the bulk, the last set of equations is not closed: 
one first needs to solve the collective dynamics for 4>q(t) 
before one can compute (j)q{t). As it should, in the 
nm ^ and rif — ^ limits, the equations for the tagged- 
particle dynamics in bulk systems [14] and those for the 
collective dynamics in the zero fluid density limit [19] arc 
recovered, respectively. The same remarks as in Ref. [19] 
apply: the presence of a linear term in the memory ker- 
nel, the disappearance of any explicit reference to the ma- 
trix, the presence of connected and disconnected struc- 
tural quantities only, and the fact that v^]^ is the same 
as in a bulk system, with connected quantities simply re- 
placing the fluid structure factor and direct correlation 
function. We refer the reader to the previous work for 
the corresponding discussion. 

Finally, from the above formulas, one can derive an 
equation for the ti me evolution of th e mean-squared dis- 
placement Sr'^{t) = — r'^(O)l^), by using the small 
wave vector behavior of the tagged-particle density cor- 
relation function, (j>q{t) = 1 — q'^Sr^ {t) /6 + 0{q'^) . It reads 
[note that it is possible to integrate once, thanks to the 

initial conditions Sr^{0) = and 6r'^{0) = 0] 

Sr2{t) + ^ f drm^^'^it - r)6r\T) = (14) 

nis Jq nis 

with mMSD(t) = liuiq^oq^mlit) = "f6{t) +m^^^^^^^\t), 
where 7 = limg^o Q^^l relates to the fast dynamical pro- 
cesses and where the mode-coupling part of the memory 
kernel is given by 
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CONNECTIONS WITH PREVIOUS 
THEORIES AND A VARIANT 



The physics of QA mixtures involves a number of 
subtleties, related, for instance, to the splitting of the 
static correlations into connected and disconnected parts. 
These might induce special difliculties and put con- 
straints on the proper way to set up a dynamical theory 
for these systems. In this technical section, we illustrate 
how they manifest themselves within the MCT frame- 
work, in an attempt to provide guidelines for future work 
based on other formalisms (see Ref. [33], for instance). 



In Ref. [19], the MCT equations for the collective dy- 
namics in a QA binary mixture have been compared to 
those for the residual dynamics of a bulk system in its 
ideal glassy phase. Strong analogies have been found in 
the overall structure of both theories, which originate 
in the presence of time-persistent density fluctuations 
in both cases. But there are significant and irreducible 
diff'erences as well, which reflect the different origins of 
these frozen fluctuations, static and disorder-induced in 
the former case, dynamical and self-induced in the latter. 

The main conclusion from this study was that, in gen- 
eral, it does not seem sensible to try and derive a dynam- 
ical theory for QA systems by simply taking the limit of 
the corresponding approach for fully annealed mixtures 
in which one component representing the solid matrix 
would become immobile, or at least extreme care should 
be taken. Since this contradicts naive expc;c;tations, we 
begin the present section with a similar analysis of the 
MCT for the tagged-particle dynamics. This is clearly 
a simpler problem, since by construction no collective 
dynamical phenomena come into play and there is no 
disconnected component in the tagged-particle density 
correlation function. 

In order to complete the discussion of Ref. [19], we 
first consider the case of a one-component bulk fluid in 
its ideal glassy phase. The tagged-particle dynamics as 
expressed in terms of the residual relaxation of the fiuid 
density fluctuations obeys mode-coupling equations of 
the same form as above, with vertices [14] 
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{(l-/|q-k|)'S'fq_k|}, 

(17a) 



iq-k|] {/|q-k|'S^-k|} 



(17b) 

where nf is the density of the fluid, its structure fac- 
tor, (?J the Fourier transform of the single-particle-fiuid 
direct correlation function, and fq the Debye- Waller fac- 
tor of the glass. To make contact with the theory for 
QA systems, we now need an expression for /if/*-^-* in 
Eq. (13b). Since the present problem corresponds to a 
situation where there are fluid density fluctuations only, 
we shall take advantage of the fact that the matrix does 
not appear explicitly in the MCT equations for QA sys- 
tems and use the Ornstein-Zernike (OZ) equations for 
a simpler example of a fiuid in a random environment, 
namely, a fiuid plunged in a Gaussian random field [34] . 
One then gets [35, 36] 



ssf(c) 



(18) 

where c^^^^ is the Fourier transform of the single-particle- 

single-particlc blocked direct correlation function. Thus, 
we observe that there is an exact correspondence between 
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the two theories, with the same role played by Sq^'^^ and If this result is confronted to the relevant expression of 



(1 — fq)Sg on the one hand, Sq^^^ and fqS^ on the other 
hand, provided one disregards in Eq. (18) the terms in- 
volving blocked direct correlation functions. This in par- 
ticular insures that the full and connected direct correla- 
tion functions coincide. 

As a second step, let us now approach the problem 
following the empirical strategy outlined above, i.e., we 
shall consider the QA binary mixture as a special limit 
of a fully annealed binary mixture. As already pointed 
out, this approach is ineffective for the description of the 
collective dynamics [19]. For a tagged particle moving in 
a fully annealed binary mixture, the MCT result for the 
memory kernel reads (for simplicity, we keep the same 
species labels as in the QA binary mixture) [22, 24, 37] 
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where c^™ is the Fourier transform of the single-particle- 
matrix direct correlation function and, in the definitions 
of the density correlation functions F^^{t), Fg'"(t), and 
F^™{t), the usual normalization for binary mixtures has 
been used, hence the prefactor rif + rim. We might now 
adapt this equation to the case of a QA binary mixture, 
for which Fg{t) splits into relaxing and frozen parts ac- 
cording to 

(nf + n^)if (t) = nf^f*^) + rif^f ^^^(i), (20a) 
{t) are actually time-independent 



while F^'^it) and F^ 
and given by 
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Sq and S™™ are the fluid-matrix and matrix-matrix 

structure factors, respectively, and, like Sq'"^^ and Sq^'^\ 
they are defined with the normalization suitable for QA 
mixtures (see Appendix A). One then obtains mode- 
coupling equations of the same form as above, with ver- 
tices 
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an exact correspondence with equations (13) is once 
again found, under the same conditions as in the pre- 
vious calculation. 

So, we find that, in both cases, there would be a per- 
fect consistency between the derived equations and the 
theory developed for fluids in random environments if 
all blocked direct correlation functions were identically 
zero. Or, stated differently, both approaches would be 
perfectly acceptable starting points for heuristic deriva- 
tions of the present theory if this condition was fulfilled. 
Obviously, it is not, and this provides further arguments 
in favor of a differentiated treatment of fully annealed 
and quenched- annealed mixtures. 

However, at the same time, the approximation which 
consists in neglecting the blocked direct correlation func- 
tions is a very common one, which for historical reasons 
is often referred to as the Madden-Glandt approximation 
[28, 29, 30]. Its widespread use originates in the fact that 
these functions are among the specific features of fluids 
in quenched disordered environments which are the most 
difficult to capture with simple approximations, and, for 
instance, many standard closures of the replica OZ equa- 
tions erroneously prescribe that they vanish identically 
[28, 29, 30, 38]. It is thus hardly surprising that they 
appear in the errors made when trying to infer a theory 
for QA systems starting from a similar approach for fully 
annealed systems, for which these functions are meaning- 
less. 

The crucial simplification in the case of the tagged- 
particle dynamics, compared to the collective dynamics, 
is that the Madden-Glandt approximation seems to be 
the only required approximation [39]. This is explicitly 
demonstrated here in the framework of the MCT, but we 
expect that analogous situations could occur with other 
theoretical schemes. So, the conclusion could be that it 
does not appear unreasonable to try and derive empirical 
theories for the tagged-particle dynamics in QA systems 
starting from approaches developed for fully equilibrated 
systems, but this will usually be at the cost of additional 
implicit approximations like Madden-Glandt's. Fortu- 
nately, most of the time, such approximations have quan- 
titative consequences only and the cost appears modest. 
Note however that one serious restriction remains: since 
no such simple empirical scheme seems to be gcncrically 
applicable to the collective dynamics [19], it will always 
be difficult to give a proper account of the situations in 
which this part of the dynamics matters, for instance, if 
the coupling of the single-particle and collective dynam- 
ics is strong. 

Another illustration of the subtle interplay between the 
dynamical theory and the peculiar structure of the static 
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correlations in fluids evolving in random environments is 
provided by the alternative derivation of the MCT for 
the tagged-particle dynamics in QA systems which has 
been discarded in the previous section. It is based on the 
observation that the only forces exerted on the tagged 
molecule are those due to the fluid and the random ma- 
trix. So, on simple physical grounds, one expects that the 
only mode-coupling contributions to the relaxation ker- 
nel will come from the products Pi^Sp^^ ^^, Pk(/'q-k)' 
PkPq-k- Projecting the random forces on this set of vari- 
ables using standard projection operator methods, one 
obtains a memory kernel which is essentially the same as 
above, except that one has to replace in Eq. (13b) ^|*q^^k| 

2 
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as given by Eq. (22) by fh^^^l^^ 

An immediate reason to reject this approach is that, 
in the limit rif — > where the dynamics should be de- 
termined by the properties of the matrix and of the 
tagged particle only, an unphysical dependence on the 
fluid would remain in the theory because of the last term 
of this expression. This is avoided with the derivation 
retained in the previous section, which is thus clearly 
preferable. 

The only difference between the two approaches is that, 
by first working with a finite concentration of tagged par- 
ticles, the correct one incorporates contributions associ- 
ated with their frozen density fluctuations, which seem 
to be preserved when the concentration is sent to zero. 
In this respect, it is pictorial that the two approaches 
happen to agree when the fluid and tagged particles are 
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identical, thanks to the equalities cl 
Indeed, it is in this unique case that the frozen fluid 
and tagged-particle density fluctuations share the same 
statistics and that the former can act in the incorrect ap- 
proach as substitutes for the latter which would otherwise 
be missing. Thus, the inclusion of contributions due to 
the frozen tagged-particle density fluctuations, which arc 
rather elusive as far as species with vanishing concentra- 
tion are concerned, appears as an important requirement 
in order to derive an acceptable theory. Retrospectively, 
failure to take into account similar contributions in spe- 
cial limiting cases (vanishing fluid density, adsorbed ideal 
gas) of the collective dynamics surely explains some am- 
biguous results pointed out in Ref. [19]. At present, the 
physics, if any, behind this constraint is unclear. 



IV. DYNAMICAL PHASE DIAGRAMS 

We now turn to the quantitative results of the theory, 
which require numerical solutions of the MCT equations. 
In the present section, we complete the dynamical phase 
diagrams of Ref. [19], where the liquid-glass transition 
lines for two simple model systems had been determined, 
while, in the next one, the full dynamical scenarios for 
the tagged-particle motion will be discussed. 

The systems to be considered arc two closely related 
QA mixture models in which both the fluid-fluid and 



fluid-matrix interactions are pure hard core repulsions of 
the same diameter d. The only difference between them 
lies in the matrix correlations. In model I, the matrix 
particles are not allowed to overlap, so that the configu- 
rational statistics of the porous medium is that of an equi- 
librium hard sphere fluid, while in model II the matrix 
particles are completely uncorrelated and overlap freely, 
as in an ideal gas. For both systems, the two dimension- 
less densities (/if = Trn[(P/6 and = Trn„^cP/6 will be 
used as control parameters, and the structural quantities 
required by the dynamical theory will be computed with 
the Percus-Ycvick (PY) approximation [28, 29, 30, 38]. 

For simplicity, we shall concentrate on the case where 
the tagged particle is one of the fluid particles, so that 
in Eqs. (13) and (16) one has to take c^q^'^'^ = Cq^'^^ and 
^ss(b) _ ^ff(b) ^ 'j'jjgjj^ ^^j^g tagged-particle dynamics coin- 
cides with the self part of the collective dynamics. 

From the point of view of the single-particle dy- 
namics, two distinct phases exist in MCT. If for all q 
limt^ao <Pq{t) = 0, one is in the diffusive regime, where 
the tagged particle can explore an inflnite domain and 
move infinitely far away from its initial position. If on 
the contrary for all q limt_+oo 4'qi't) = fq 0, one is in 
the localized regime, where the tagged particle is trapped 
in a finite domain. 

The single-particle nonergodicity parameter /|, which 
is usually called the Lamb-Mossbauer factor, is the solu- 
tion of the nonlinear set of equations 



J a 



(27r)3 



(2) 



/fe/|c 



(23) 



where fg is the collective dynamics nonergodicity param- 
eter or Debye- Waller factor computed in Ref. [19]. This 
equation has to be solved numerically in order to locate 
the diffusive and localized phases when (j)f and 0m are 
varied. All computations in the present paper have been 
achieved using the methods of Refs. [40, 41], to which the 
interested reader is referred for technical details, with the 
same parameters as in Ref. [19]. As pointed out in the 
latter work, the discretization used for the wave vector 
integration in Eq. (23) involves a cutoff of the low q diver- 
gence of the memory kernels. This has the advantage of 
eliminating spurious long time anomalies originating in 
this divergence [42, 43, 44], but at the cost of quantitative 
uncertainties due to the arbitrariness of the cutoff value. 
They have been estimated in two ways, first by test cal- 
culations on a finer q grid, second by a comparison with 
the results obtained within an additional hydrodynamic 
approximation (see Appendix C) which allows one to in- 
tegrate exactly over the full q range [26]. It has been 
found that the errors on the location of the transition 
points never exceed a few percents. 

The dynamical phase diagrams of models I and II are 
reported in Fig. 2. The liquid-glass transition lines of 
types A and B are those obtained in Ref. [19]. As dis- 
cussed there, for model I they join smoothly at a degen- 
erate higher-order A3 singularity, while for model II a 
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FIG. 2: Dynamical phase diagrams of two types of hard-sphere quenched-annealed binary mixtures, (a) Model I: non- 
overlapping matrix particles, (b) Model II: freely overlapping matrix particles. For the collective dynamics, the A3 endpoint 
singularities are denoted by E and the crossing poiut between the type A and B transition lines in model II is denoted by C. 
For the tagged-particle dynamics, the crossing points between the diffusion-localization and glass transition lines are denoted 
by T (for triple point). 



crossing point followed by a line of glass-glass transitions 

terminating at an ordinary ^3 singularity is obtained. 
Except for these details, both phase diagrams are very 
similar at the level of the collective dynamics. 

This remains true when the tagged-particle dynamics 
is considered. For both models, two types of transitions 
between the diffusive and localized phases are found. 

At low matrix densities, the glass and localization tran- 
sitions occur simultaneously. In fact, it is the spon- 
taneous arrest of the collective dynamics which trig- 
gers the arrest of the tagged-particle dynamics through 
their coupling in the bilinear term of the memory ker- 
nel (12). This is the well known scenario found in all 
MCT studies of one-component bulk systems so far, in 
which the tagged-particle dynamics inherits its discontin- 
uous (type B) character and quantitative features (criti- 
cal laws, exponent parameter A) from the collective dy- 
namics [14, 41]. Its extension to weakly confined systems 
is easily understood from simple continuity arguments in 
the vicinity of the bulk transition at 0i„ = 0. 

These arguments break down at higher matrix densi- 
ties, where one finds that the tagged particles are already 
localized when the ideal glass transition occurs. Thus, 
the latter only manifests itself as a secondary dynamical 
transition in an already arrested single-particle dynam- 
ics, with characteristic features imposed by the collective 
dynamics just like in the previous scenario. In fact, in 
this density domain, the liquid-glass transition is always 
preceded by a continuous diffusion-localization transition 
driven by the linear term of the memory kernel (12). 
These transitions, which only manifest themselves in the 
single-particle dynamics, are special cases of type A tran- 



sitions, characterized by an exponent parameter A = 0, 
like the pure diffusion-localization transition at 0f = 
[19, 25, 26]. Similar ones have been previously discussed 
in the framework of schematic models [14, 45, 46, 47], 
and, more interestingly, have also been found in early 
MCT studies of bulk binary hard sphere mixtures with 
large size asymmetry, where they take place in the ideal 
glassy phase and correspond to the localization of the 
smaller species in the voids of the glass [22, 23, 24]. Su- 
perficially, the analogy between the two systems can be 
understood by seeing the disordered porous solid in the 
QA mixture as a glassy component or by considering the 
induction of permanent fluid density fluctuations by the 
quenched random environment as some kind of previtri- 
fication. But we stress once again that, as illustrated in 
Ref. [19] and in the previous section, such simple pictures 
do not in general lend themselves to the development of 
rigorous approaches to the dynamics of QA systems. 

As seen in Fig. 2, the corresponding transition line 
starts on the <p[ = axis at the same point as the liquid- 
glass transition line, since in this limit the collective dy- 
namics exactly reduces to the single-particle dynamics. 
As is increased, it runs across the fluid domain as ob- 
tained from the study of the collective dynamics, and 
it terminates when it intercepts the type B liquid-glass 
transition line at some finite value of <^in. 

The final result is thus that the dynamical phase dia- 
grams of both Q A mixtme models display three domains, 
as shown in Fig. 2. For low overall densities, the dy- 
namics is fluid-like, with both the collective and tagged- 
particle density fluctuations relaxing to their equilibrium 
values at long times. For intermediate densities, the dy- 
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namics is localized, with a complete relaxation of the col- 
lective density fluctuations only. And for high densities, 
glassy dynamics is found, in which neither the collective 
nor the single-particle density fluctuations do return to 
equilibrium. These three domains meet at two points: 
first, at the pure diffusion-localization transition point at 

= 0, because, as mentioned above, there coincide the 
collective and tagged-particle dynamics; second, at the 
crossing point between the liquid-glass and continuous 
diffusion-localization transition lines. For definiteness, 
because of the three phases, we shall call this point a 
triple point denoted by T in Fig. 2. But it should be 
remembered that phase coexistence is excluded in MCT, 
because of a maximum property of the physical solution 
of the MCT equations [14]. It results that point T cor- 
responds unambiguously to a glassy state. 

The clear separation between the diffusion-localization 
and liquid-glass transition lines allows one to completely 
disregard the collective dynamics when dealing with the 
asymptotic properties of the tagged-particle dynamics 
in the vicinity of the former transition. This possibil- 
ity is particularly interesting if an additional hydrody- 
namic approximation is performed, since significant ana- 
lytic progress can then be made in the study of the con- 
tinuous diffusion-localization transition. This approach, 
which extends Leutheusser's theory [26] to finite fluid 
densities, is discussed in Appendix C. 

Beside the nature of the domains and transition lines 
met in the phase diagrams, their shape is of interest as 
well. A remarkable result of Ref. [19] has been the pre- 
diction of a reentry phenomenon in the liquid-glass tran- 
sition line for high matrix- low fluid densities, which has 
been interpreted as the signature of a dynamical decorre- 
lation mechanism due to fluid-fluid collisions c;ompeting 
with the reduction of free volume and the localization 
effect due to the disordered porous solid. A similar fea- 
ture, which might be understood using the same physical 
arguments, is also present in the continuous diffusion- 
localization transition line, but with a much weaker am- 
plitude, so that one needs the strong magnifications of 
Fig. 3 in order to visualize it. 

As already pointed out in Ref. [19], a word of caution 
is needed about this finding. Indeed, for systems ruled 
by hard core fiuid-matrix interactions like those studied 
in the present work, a number of recent extensive sim- 
ulation studies [48, 49, 50, 51, 52, 53, 54] have clearly 
demonstrated that the localization transition is actually 
driven by the percolation transition of the matrix void 
space, i.e., localization occurs because, above a certain 
critical matrix density, the void space only consists of 
finite disconnected domains. Such a scenario definitely 
rules out the possibility of a reentry phenomenon in the 
diffusion-localization transition line as obtained above. 
Indeed, if at a given value of 0m a system with (pf = 
is localized because of the onset of the percolation tran- 
sition, any system with a finite (pf at the same matrix 
density will be localized as well, since a variation of (j)f 
obviously has no effect on the topology of the matrix 
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FIG. 3: Magnified views of the continuous diffusion- 
localization transition lines for high matrix densities. A reen- 
trant behavior is clearly visible, (a) Model I, (b) Model II. 



empty space. 

Fortunately, the amplitude of the reentry phenomenon 
predicted by the MCT for the diffusion-localization tran- 
sition is very small, so the inconsistency is quantitatively 
not too serious, but, clearly, this prediction should not 
be taken literally. Instead, a reasonable expectation is 
typically a nonmonotonic variation of the self-diffusion 
coefficient in the slowly diffusive regime just below the 
localization threshold, with first an increase due to fluid- 
fluid collisions when increasing ([>[ from zero. In support 
of this suggestion, we recall that, in a number of systems 
for which the MCT has predicted reentrant transition 
lines, such a reentry has been observed in the isodiffusiv- 
ity curves as obtained by molecular dynamics simulations 
[55, 56, 57, 58, 59, 60, 61]. The same behavior should also 
be expected in systems with soft repulsive fluid-matrix 
interactions, for which the percolation concepts are only 
approximate. In this respect, it is encouraging that this 
is precisely what has been observed in a computer simu- 
lation study of a two-dimensional lattice gas model with 
fixed randomly placed penetrable scatterers [62]. 



V. DYNAMICAL SCENARIOS 

According to the phase diagrams reported in the pre- 
vious section, both QA mixtures considered in this work 
display three generic types of global dynamical scenarios 
in response to increases of the overall density of the fluid- 
matrix system. For small (j)^, a bulk- like type B ideal 
glass transition scenario prevails, in which both the col- 
lective and single-particle density fluctuations simultane- 
ously show discontinuities. For moderate (j)^, a continu- 
ous diffusion-localization transition first occurs, followed 
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at a higher density by a type B ideal glass transition. 
Finally, for high (f>^, the latter transition is of type A. 

It is the purpose of this section to provide illustra- 
tions of these three typical scenarios and to discuss their 
features which might be relevant for comparisons of ex- 
perimental or simulation data with the predictions of the 
theory. As in Ref. [19], the analytic asymptotic results 
shall be quoted without their proofs, which can be found 
in or adapted from Rcf. [14] (sec also Rcf. [47] for type A 
transitions). Note that these results, which are obtained 
under the assumption of regular vertices, as well as the 
method of solution of the mode-coupling equations used 
in this work [40, 41], which involves a cutoff of the mem- 
ory kernels at low q, cannot give an account of possible 
dynamical features originating in the divergence of the 
mode-coupling vertices when g — » [42, 43, 44]. While 
this might look unsatisfactory from a mathematical point 
of view, such an approximation is in fact required for 
physical reasons, since the low q singularity is actually 
known to be an ill feature of the mode-coupling approx- 
imation which results in spurious long time anomalies 
[43]. 

For simplicity and consistency with Rcf. [19], wc fol- 
low Refs. [40] and [41] and consider the dynamics in the 
overdamped limit valid for Brownian systems. The gen- 
eralized Langevin equations (11) and (14) then reduce 
to 

r'M*) + '^^(i) + /* S^^''^ it - ^)<^g(^) = (24) 
Jo 

and 

Sr'' (i) + Do f dT m^SD(MC) _ ^-,^^2 ^ g^,^^^ ^25) 
Jo 

with Do the short-time diffusivity, t| = l/{Doq^), and 
the initial conditions </)^(0) = 1 and (5r^(0) = 0. This 
simplification changes the short time transient part of 
the dynamics, but not its long time properties. In the 
following, the particle diameter d is used as the unit of 
length and the unit of time is chosen such that Do = 
1/160. Also, since both QA models are found essentially 
equivalent, only model I will be considered. 

For the bulk-like scenario found for small 4>m, it is now 
well known from extensive studies of one-component bulk 
systems [14, 41] that, in this case, the essential features 
of the tagged-particle dynamics are simple reflections of 
those of the collective dynamics, since it is the latter 
which actually drives the system through its unique er- 
godicity breaking transition. So we do not repeat the de- 
scription of the type B scenario given in Ref. [19], which 
remains valid in the present case without change, and we 
omit the analogues for /| and (f)q{t) of Figs. 5(a), 7(a), 
and 8(a) of Ref. [19], in which the evolutions with density 
of fq and 0q(t) were reported, since they would display 
exactly the same patterns. Instead, we concentrate on 
the distinctive aspects of the single-particle dynamics, 
i.e., the wave vector dependence of the Lamb-Mossbauer 
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FIG. 4: Wave vector dcpcii(iciice of the single-particle 
nonergodicity parameter or Lamb-Mossbauer factor /| = 
limt^oo <t>qit) for model I at matrix density 0m = 0.05. Prom 
left to right, bottom to top: 4>f = (jff, 1.01(/l)f, l.O50f , and 
1.10^ 

factor and the time dependence of the mean-squared dis- 
placement. As in Ref. [19], we choose (j)^ = 0.05 for the 
present illustrations. 

Fig. 4 shows as a function of q at different fluid 
densities in the glassy phase. As in bulk systems, the 
curves start at 1 for q — Q, arc approximately Gaussian- 
shaped and widen rapidly when 0r is increased. For a 
given distance to the transition, they are generically nar- 
rower than the corresponding curves for the bulk, reflect- 
ing the fact that, when increases, the system evolves 
towards a continuous glass transition scenario. 

The mean-squared displacements at different fluid den- 
sities near the critical one arc reported in Fig. 5. For 
small t, short-time diffusion with diffusivity I^o is ob- 
served. Then, in the ergodic phase, the familiar two-step 
dynamics typical of the discontinuous ideal glass tran- 
sition scenario sets in, with, as 0f is approached, the 
development of a plateau reflecting transient localiza- 
tion, which crosses over for large t to the diffusion be- 
havior 6r^{t) = 6Dt. The long-time diffusion coefficient 
is generically given by [41, 63] 

D = joo^ . (26) 

1 + Do m^s°(^c)(t)dt 
Jo 

and thus vanishes in the present scenario according to 
the a relaxation scaling law, i.e., cx |(/)f — 0^1^/^"+^/^'', 
where a and b represent the critical decay and von Schwei- 
dler exponents of the collective dynamics, respectively. 
When the transition density is reached and exceeded, 
the plateau lasts forever, the tagged particle being per- 
manently trapped. Close to the bifurcation point, the 
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FIG. 5: Time evolution of the mean-squared displacement 
Sr^{t) for model I at matrix density 0m = 0.05, in the vicin- 
ity of the liquid-glass transition threshold (f>(. Prom left to 

right, top to bottom: (^f = 0.9<j!>f, 0.99<?!-f , 0.999<^f, 0.9999(|!>f , 
0.99999<^f, 1.00001(^f, 1.0001<^f, l.OOl^f, l.Ol^f, l.l(j>l The 
horizontal dashed line marks 6rf^, with rf the critical local- 
ization length. The dotted line shows the short-time diffusive 
behavior 5r'^{t) = GDot. 



FIG. 6: Time evolution of the single-particle density corre- 
lation function (j>\{t) at g ~ 7.09/d for model I at matrix 
density (/>m = 0.1, in the vicinity of the localization threshold 
(f)\. From left to right, bottom to top: (j>t — 0.94>\, 0.990', 
O.9990f, 0.99990^, 0.999990,, l.OOOO10f, l.OOO10f, l.OO10f, 
l.O10f, 1.10f. The dashed line shows the analytically derived 
critical decay law. 



localization length r\, defined through 



so that 



lim Sr^t) = 6rf, 

t—*-oo 



/■ + ~ k^dk ( (2),,, ^ 



(27) 



(28) 



displays a square-root singularity, r\—ry oc — (0>j— 0^)'^/^, 
with the finite loealization length at the transition. In 
both the diffusive and localized phases, the dynamics in 
the vicinity of the plateau is described by the same scal- 
ing functions as the density fluctuations, with a critical 
amplitude oc |(/)f — 1^^^ and a characteristic time scale 

The two other scenarios, occiiring for moderate and 
high ^m, have in common to start with a continuous 
diffusion-localization transition which affects only the 
tagged-particle dynamics and leaves the collective quan- 
tities free of any singularity. It will now be discussed in 
some detail and illustrated with computations at = 
0.1. 

The time evolution of the single-particle density corre- 
lation function 0q(t) at g ~ 7.09/d, corresponding to the 

main peak of S^^'^^ , is reported in Fig. 6 in a log- log plot 
(the curves for other values of q are qualitatively simi- 
lar). In both the diffusive and localized phases near the 
transition, a single step relaxation is found. When the 
localization threshold 0f is approached from below, the 



slowing-down of the dynamics manifests itself through 
the development of a weak long time tail which extends 
to longer and longer times as the density is increased. At 
the transition, it lasts indefinitely, and, above the thresh- 
old, it progressively recedes and continuously turns into a 
finite Lamb-Mossbauer factor. From a quantitative point 
of view, the MCT diffusion-localization transition is most 
easily understood as a type A bifurcation characterized 
by an exponent parameter A = and thus a critical decay 
exponent a = 1/2. It results that, close to the transition 
in the localized phase, /| oc ((^f — <^{:). A reduction theo- 
rem also holds, stating that, in both phases and for long 
enough times, the wave vector and time dependence of 
(lfq{t) factorize according to 



\{t) - lim 4.\{t) = h\G\t). 



(29) 



At the threshold, one gets the critical decay 

G%t) = {tl/tf'\ (30) 

where is a time scale obtained by matching the short 
and long time dynamics, while, for finite values of |0f — 
4>\\, B. scaling law is found, 



(31) 



with oc — and /t^ oc — g^{f) is 

discussed in detail in Ref. [14] and one finds in particular 
that 



g%t) = t 



- +-1/2 



(32) 
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FIG. 7: Time evolution of the mean-squared displacement 
(5r^(t) for model I at matrix density <^m = 0.1, in the vicin- 
ity of the localization threshold From left to right, top 
to bottom: 4>t = 0.9(/>f, 0.m<t)\, 0.999<^f, O.99990f, 0.99999(?if, 
1.00001(^f, l.OOOl^f, 1.001(?if, l.Ol^f, l.l<^f. The dotted and 
dashed lines show, respectively, the short-time diffusive be- 
havior 5r^{t) = QDot and the analytically derived critical 
subdiffusive law 5r'^{t) oc t^^^. 

for small t, so that the critical decay law is reached con- 
timioiisly when (j)[ (p\- Note that the sign of 4>{ — <p\ 
is irrelevant in the dynamical scaling laws. An impor- 
tant feature of all type A scenarios follows, which is the 
symmetric departure from the critical decay law (30) at 
long times for state points located in the ergodic and 
nonergodic phases at the same distance from the transi- 
tion. All these analytic predictions are easily confirmed 
by inspection of Fig. 6. 

Since, according to Eqs. (14)-(16), the time evolution 
of the mean-squared displacement is controled by that 
of the single-particle density fluctuations, the anoma- 
lous dynamics of the latter near the diffusion-localization 
transition is naturally transfered to the former. Indeed, 
as seen in Fig. 7, the long time tail of (jfq{t) induces 
the development in 6r'^{t) of a subdiffusive regime at 
intermediate times, where dr^^t) oc t^/^. This regime 
would persist forever exactly at the transition. For fi- 
nite (f)f — it ends at a time of the order of r^, where 
it crosses over to ordinary diffusion 5r'^{t) = 6Dt in the 
diffusive phase or to a constant Sr'^{t) ~ 6rf in the local- 
ized phase. When the transition density is approached 
from below, the long-time diffusion coefficient vanishes 
linearly, D oc \(l)f — ^f\, while the localization length di- 
verges according to r\ oc |0f — 0^ j^-'^/^ when the threshold 
is approached from above. Here again, all these asymp- 
totic laws can be easily checked in the figure. 

Once localization has occured, further increase of the 
overall density drives the system towards its ideal liquid- 



glass transition, which can be either discontinuous or con- 
tinuous. The first case is realized for 0^ = 0.1, which was 
considered above from the point of view of the diffusion- 
localization transition. To illustrate the second case, we 
will choose (j)f = 0.15 and use (fim as the external control 
parameter, since, with the QA mixture models studied in 
this work, it is not possible to find a suitable constant (j)ni 
trajectory in the dynamical phase diagrams. Note that 
this exchange of the roles played by c/if and (/),„ does not 
affect the general critical properties of the theory, pro- 
vided they are formulated in terms of the relevant control 
parameter. These critical properties will not be discussed 
in detail. Indeed, it can be very easily demonstrated, by 
introducing shifted and rescaled tagged-particle density 
correlation functions so as to eliminate the contributions 
due to localization, that there is actually no formal differ- 
ence between the present situation and the simpler case 
of the bulk-like scenario. So, here as well, the anoma- 
lous dynamics of the tagged-particle density fluctuations 
in the vicinity of the ideal liquid-glass transition is com- 
pletely controled by the collective dynamics, which im- 
poses the main qualitative and quantitative features of 
the relaxation. 

The complete double transition scenario at 0m = 0.1 
is illustrated in Figs. 8, 9, and 10. The data correspond- 
ing to the lowest densities provide further illustrations 
of the features mentioned above for the dynamics near 
the diffusion-localization transition threshold (p\. This 
includes the continuous development of the long time tail 
and finite asymptote in (tfq{t), now shown in a semi-log 
plot in Fig. 8, and the continuous and linear increase of 
/| from zero at i/ij-, as seen in Fig. 9. Then, approaching 
the glass transition density the single-particle den- 
sity correlation functions start to display a two step de- 
cay, which is clearly visible in Fig. 8 and is naturally ex- 
pected from a type B glass transition scenario. When the 
second relaxation step disappears because of the spon- 
taneous arrest of the collective dynamics, a discontinu- 
ity followed by a square root singularity occurs in the 
Lamb-Mossbauer factor, as seen in Fig. 9, whose inset 
shows that as a function of q always remains bell- 
shaped. All these anomalous behaviors are transmitted 
to the mean-squared displacement through the memory 
kernel m^^^^'^^''^\t) and result in the two step relaxation 
pattern visible in Fig. 10. 

The analogous data for the double transition scenario 
at 0f — 0.15 are reported in Figs. 11, 12, and 13. Com- 
pared to the previous case, there is no significant differ- 
ence from the point of view of the dynamics near the 
diffusion- localization transition threshold But, as 

the matrix density is increased, the relaxation of the 
tagged-particle density fluctuations always remains sin- 
gle stepped, as shown by Fig. 11. The type A liquid-glass 
transition signals itself by a considerable slowing down of 
the approach by (jf^ {t) of its infinite time value, which pro- 
cedes via an algebraic decay law with the same exponent 
as the collective dynamics. The single-particle nonergod- 
icity parameter /| remains continuous at the transition. 
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FIG. 8: Time evolution of the single-particle density correla- 
tion function at g ~ 7.09/d for model I at matrix den- 
sity (pm = 0.1, in the vicinity of the diffusion- localization and 
liquid-glass transition thresholds (f>\ and (f>f. From left to right, 
bottom to top: <f>{ = 0.9(^f, 0.99(j>\, 1.01(j)\, ~ O.90f, 

0.99<^f, 0.999<^f, O.99990f, 0.99999<?!)f , 1.00001<?!)f , l.OOOl^f, 
l.OO10f , l.O10f , 1.10J. The lower and higher horizontal 
dashed lines mark the Lamb-Mossbauer factors in the lim- 
its 0f — »■ (j>f~ and 0f —> ct>f~^ , respectively. 
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FIG. 10: Time evolution of the mean-squared displacement 
(5r^(t) for model I at matrix density 0m = 0.1, in the 
vicinity of the diffusion-localization and liquid-glass transi- 
tion thresholds (^J and From left to right, top to bot- 
tom: 4>i = 0.9(^f, 0.99(^f, 0.99901, O.99990f, 1.0001^, l.OOl^f, 
1.014, 1.10} ~ O.90f, O.990f, O.9990f, O.99990f, O.999990f , 
l.OOOO10f , l.OOO10f , 1.00101, l.O10f , 1.10f . The higher and 
lower horizontal dashed lines mark 6rf in the limits 0f ^ 0f ~ 
and 0f 0f^, respectively. The dotted line shows the short- 
time diffusive behavior (t) = ^D^t. 
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FIG. 9: Density dependence of the single-particle nonergod- 
icity parameter or Lamb-Mossbaucr factor /| = limt^oo 0q(i) 
at g ~ 7.09/(i for model I at matrix density 0m = 0.1. Vertical 
dotted lines and diamonds indicate the diffusion-localization 
and liquid-glass transitions. Inset: Wave vector dependence 
of /°. From left to right, bottom to top: 0f = l.O50f, 1.10f, 
0f", 0f+, l.O50f, and 1.10f. 



which is only marked by an abrupt but finite change of 
slope visible in Fig. 12. Finally, here as well, the time evo- 
lution of the mean-squared displacement br^ (t) , reported 
in Fig. 13, is directly influenced by this anomalous relax- 
ation of the density fluctuations and, as a result, shows 
a regime of continuous critical dynamics in the localized 
phase near the glass transition point. 

Like in Ref. [19], we shall leave aside the peculiar dy- 
namics found in the vicinity of the state points corre- 
sponding to higher-order singularities or crossing points, 
either fully embedded in the liquid-glass transition lines 
(points E and C) or at their junctions with the diffusion- 
localization transition lines (points T). They certainly 
are interesting, but such a study would be quite techni- 
cal, requiring fine-tuning of the parameters of the models 
and refined mathematical tools [14, 64, 65, 66, 67]. So 
we prefer to rather concentrate on the generic scenarios, 
which are the most likely to be observed in experiments 
or computer simulations, and wc keep such an advanced 
discussion for future work, when it will be useful, for in- 
stance, if relevant data become available. 



VI. DISCUSSION AND CONCLUSION 

In this paper, a MCT for the slow dynamics of a tagged 
particle moving in a fiuid adsorbed in a disordered porous 
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FIG. 11: Time evolution of the single-particle density cor- 
relation function 4>q{t) at g ~ 7.09/d for model I at fluid 
density <f>f = 0.15, in the vicinity of the diffusion-localization 

and liquid-glass transition thresholds <f)\^ and 01^. From left 
to right, bottom to top: (pm = 0.901^, O.990|n, LOl^l^i, 1.10m, 
0.901,, 0.9901,, 0.99901,, 1.00101,, 1.0101,, 1.101,. The hori- 
zontal dashed line marks the Lamb-Mossbauer factor at 0|,. 
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FIG. 13: Time evolution of the mean-squared displacement 
Sr^{t) for model I at fluid density 0f = 0.15, in the vicin- 
ity of the diffusion-localization and liquid-glass transition 
thresholds 0^ and 0|,. From left to right, top to bottom: 
0„, = O.90L,, O.990|„, O.9990l„, O.99990|„, l.OOO10|„, 1.0010^, 
1.O10L,, 1.10L,, 0.901,, 0.9901,, 0.99901,, 1.00101,, 1.0101,, 
1.101,. The horizontal dashed line marks 6rf at 0|,. The dot- 
ted line shows the short-time diffusive behavior 5r^{i) = 6Dot. 
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FIG. 12: Density dependence of the singlc-particlo nonergod- 
icity parameter or Lamb-Mossbauer factor fg = limt^oo 
at g ~ 7.09/rf for model I at fluid density 0f = 0.15. Vertical 
dotted lines and diamonds indicate the diffusion-localization 
and liquid-glass transitions. Inset: Wave vector dependence 
of fq. From left to right, bottom to top: 0m = 1.050m, l-10m: 
01,, 1.0501,, and 1.101,. 



solid has been developed. It complements previous work 
on the collective dynamics reported in Refs. [17, 18, 19], 
so that one might now consider as complete the mode- 
coupling theoretical framework for the study of the dy- 
namics of density fluctuations in QA systems. 

From a formal point of view, the newly derived equa- 
tions keep the features of those of Ref. [19] which were 
found appealing. Thoy arc universality, in the sense that 
the eventual dynamical equations do not contain any ex- 
plicit reference to the precise nature of the random envi- 
ronment in which the fluid evolves, and a very close simi- 
larity with other mode-coupling equations previously de- 
rived in comparable contexts. So, the overall structure of 
the theory appears quite satisfactory. Note however that 
this turns out not to be a completely trivial task to ob- 
tain such equations which display these features and also 
avoid certain subtle inconsistencies. Indeed, as shown in 
Sec. Ill, special care seems to be generically needed in 
order to properly handle the peculiar c;orrclations in QA 
systems, for reasons which are not fully clear at the mo- 
ment. In view of these difficulties, it is hoped that the 
present approach can be useful in the future as a guide 
for developments based on other formalisms. 

Prom a physical point of view, the main result of the 
present work is the prediction of a continuous diffusion- 
localization transition which can occur before and inde- 
pendently of the liquid-glass transition in situations of 
strong confinement. It is characterized by the emergence 
of various dynamical anomalies, in particular, by the de- 
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velopment of an algebraic long time tail in the tagged- 
particle density correlators and by the opening of a sub- 
diffusive regime in the mean-squared displacement. 

Dealing with the problem of a tagged particle mov- 
ing in a random medium, the possibility of a diffusion- 
localization transition is a natural expectation. However, 
the prediction that the single-particle dynamics can be- 
come singular and nonergodic while the collective dy- 
namics remains regular and ergodic might look suspect, 
owing to the fact that the former contributes to the lat- 
ter, of which it represents the so-called self part. In fact, 
the difficulty is only superficial. Firstly, it should be re- 
membered that the theory for the collective dynamics 
and the definition of the ideal liquid-glass transition are 
formulated in terms of the connected fluid density cor- 
relation functions, while the single-particle density cor- 
relation function is the self part of the total fluid den- 
sity correlation function. Because the random external 
field provided by the fixed matrix induces nonzero aver- 
age density fluctuations in the fluid at equilibrium, the 
latter, defined as 



(33) 



with S? = {pl.plJ/Ni = Sg^"^ + Sq^^\ can be written 



as 



S, 



At) 



s 



ff(b) 



(34) 



from which it is clear that (f)^{t) never decays to zero in 
a QA system. So, one might say that the fluid density 
fluctuations are never really ergodic in a QA mixture, in 
the loose sense that they saturate at a finite value, like 
in a bulk glassy system. Secondly, based on the results 
of the MCT for one-component bulk systems [14] and on 
the prevalence of rudimentary approximations like Vine- 
yard's [68], it is often believed that the collective and 
self density fluctuations should generically behave more 
or less the same way. This is not correct. In fact, one 
can easily demonstrate, using the memory function for- 
malism, that the total and single-particle density correla- 
tion functions are completely independent, formally de- 
coupled objects [69], so that the general case is rather 
the contrary, i.e., the two functions display quite dif- 
ferent behaviors [70]. Thus, despite appearances, the 
pure diffusion-localization transition scenario predicted 
in this work, which is analogous to the one predicted 
years ago for bulk binary mixtures with large size asym- 
metry [22, 23, 24], is not marred by any inconsistency. 

An interesting aspect of this prediction, that we ex- 
pect to be quite general, is that it demonstrates the 
possibility, when strong static fluid correlations [encoded 
here in /i""('^)(r)] are present in confinement, that the 
single-particle dynamics develops slow relaxational fea- 
tures which have no relation whatsoever with glassy dy- 
namics, if the latter is to be understood in the usual way 



as a collective dynamical phenomenon. We believe that 
this finding could be important for the design and the 
interpretation of experiments and computer simulations 
on the dynamics of confined glassforming systems. In- 
deed, the idea that confinement can be responsible for 
the appearence of additional slow dynamical processes 
which blend with the glassy dynamics and might obscure 
its features is definitely not novel [7, 8] (see Refs. [71, 72] 
for a specific experimental observation of such a process). 
But it does not seem to have been realized previously that 
there could be such a major difference in the way some 
of these processes are reflected in collective and single- 
particle quantities, being essentially invisible in the for- 
mer and possibly dominant in the latter. This suggests 
that, if one is mainly interested in understanding the ef- 
fect of confinement on the glass transition, it should be 
easier to concentrate on collective dynamical quantities, 
and that, if single-particle data are to be considered, one 
should be prepared to face more complex dynamical sce- 
narios which might result from the interplay between the 
glassy dynamics and confinement-induced single-particle 
relaxation processes. 

Finally, all the predictions made in the present and 
previous papers [17, 18, 19] can in principle be tested 
by computer simulations, in order to judge the quality 
of the proposed theory. As far as the single-particle dy- 
namics is concerned, this should not be too challenging, 
since, in the past few years, a significant number of simu- 
lation studies have actually been already reported, which 
precisely dealt with the tagged-particle dynamics in QA 
systems [33, 53, 73, 74, 75, 76, 77, 78, 79, 80, 81]. Un- 
fortunately, in all these studies, either the model was 
simple, but the focus was not on slow dynamics, or the 
glassy dynamics was explored, but the fluid model was 
quite complex, so that direct comparisons with the re- 
sults of the present work are never possible. We hope 
that the situation can improve very soon. 



APPENDIX A: REPLICA ORNSTEIN-ZERNIKE 
EQUATIONS 

For reference, we qiiote in this appendix the replica 
Ornstein-Zernike equations relating the various pair cor- 
relation functions mentioned in the main text. 

For the fluid-matrix part, they read [28, 29, 30, 31] 



Lmm _ -mm , -mmtm 



(Ala) 

= + n^c^h^"^ + Hfcf'^^h'^, (Alb) 

/^ff(b) ^ ^ff(b) ^ „^.fm/^fm ^ „^5ff(c)/,ff(b) + „,gff(b)/,ff(c)^ 

(Ale) 
(Aid) 



with Cq = Cq^'^^ + Cq^^^ and hq = hq^'^^ + hq^^^ . As usual, 
h and c denote total and direct correlation functions, 
respectively. / denotes the Fourier transform of / and 
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the superscripts have the same meaning as in the main 
text. 

They can be easily combined with the definitions of 
the structure factors, 



1 + nm/i™, 



(A2a) 
(A2b) 
(A2c) 
(A2d) 



(remember that = Sq'^"' + Sq'-'"' , hence = 1 + 

mh^g) to express them in terms of the direct correlation 
functions (see Appendix A of Ref. [19]). 

For the single-particle part, the relevant equations, 
which are limiting cases of those reported in Ref. [32], 
read 



^q 



Cg rig 



Hfc; 



fm 



(A3a) 



/,sf(b) = gsf(b) + ri^^f^hl^ + nfcf'^)/if ^) + nfcf 

(A3b) 
(A3c) 

(A3d) 
(A3e) 



^ss(b) 



fsf(c) _ -;Sf(c) I „,ASf(c)ffr(c) 

"-q ~ g ' ' (? q ' 



^ss(c) ^gss(c)^^^gsf(c)^sf(c)^ 



total fluid particle number is N[ = X]a=i and the 
fluid number fractions are defined as Xa = Na/N{. The 
system is at temperature T and has volume V, hence the 
total fluid density is Uf = Nf/V. Note that it is not neces- 
sary to provide any information on the disordered matrix, 
since only quantities characterizing the fluid component 
of the QA system will appear in the final results. 

The dynamical variables of interest are the time- 
dependent density fluctuations of the different fluid 
species, 



(Bl) 



where r" (t) is the position of the fluid particle j of type 

a at time t. From these quantities, one can form the 
relevant static structure factors, connected, 



(B2) 



with Sp"{t) = Pa{t) — (Pq), and blocked or disconnected. 



(B3) 



as well as the nonnormalized time-dependent connected 
density fluctuation autocorrelation functions 



with = cf^'+cf^', 
cf^'^'and = 



rsf(c) 



sf(b) 



s(c) 



APPENDIX B: MODE-COUPLING EQUATIONS 
FOR A MIXTURE ADSORBED IN A 
DISORDERED POROUS SOLID 

In this appendix, we report the MCT equations for the 
dynamics of a mixture adsorbed in a disordered porous 
matrix. They are derived following an obvious generaliza- 
tion of the procedure outlined in appendix B of Ref. [19] 
and are thus given without proof. 

The motivation for this appendix is twofold. First, 
these equations are needed in order to derive those for 
the tagged-particle dynamics, which are the main result 
of the present work. Second, most simple glassformer 
models which are considered in the literature for compar- 
isons between theoretical predictions and intensive com- 
puter simulation results are binary mixtures, for the rea- 
son that they are much less prone to crystallisation than 
one-component systems. So it is likely that sooner or 
later a version of the theory dealing with adsorbed mix- 
tures will become useful, for instance, in order to perform 
quantitative tests like the ones reported in Refs. [82, 83] 
for bulk systems. 

We consider a quenched-annealed multicomponent 
mixture containing n fluid components. Each of them 
consists of Na particles of mass nia (1 < a < n). The 



(B4) 



These functions obey standard generalized Langevin 
equations which read in matrix form, 



Fg{t)+n'gF,{t)+ dTM,(i-T)F,(r) = 0, (B5) 
Jo 

with initial conditions Fg(0) = Sq'^ and Fq(0) = 0, and 
a frequency matrix given by 



where S^'^'' ^ is the matrix inverse of S^'^-'. 

The MCT provides an expression for the kernel Mg(t), 

which reads Mq{t) = Tq6{t) + M^^'~^\t), where F, is a 
matrix of friction coefficients associated with fast dynam- 
ical processes and 



m 



aa'(2) 



(t)+m««'W(d, (B7) 



with 



(B8) 
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and 



with the generahzed diffusion coefficient Dq{z) given by 



E 

/3,/3'=l 



q,k q,k 



q-k| 



(C2) 



The vertices are given by 



c"'^^'^ 



+ 



q - (q 



a/3(l) _qjc.a/3(c) , Q " (Q - k) 1 



•|q-k|' 



q,k 



where the c, 



.a/3(c) 



+ 



rifXo 



(B9) 

(BIO) 
(Bll) 



The hydrodynamic approximation simply amounts to ne- 
glecting the q dependence of Dg{z) in Eq. (CI), where it 
is replaced by its — > limit 



D{z) 



-kBT/m.s 



z + fh^^^{z)kBT/ms 



(C3) 



s arc the Fourier transforms of the con- 



nected direct correlation functions [32]. 

The equations for the tagged-particle dynamics are 
obtained by considering a species a in the limit of 
a vanishing number fraction. One then has ^g(t) = 
Y\mx^^oF^''{t)/x„. After tedious but straightforward 
algebra, one finds the usual Langevin equation (11) and 



^s(MC)(^^ = mf\t) + mf\t), (B12) 



with 



q (q-k) 



s/3(c) .s/3'(c) ■ 



-|q-kr 



and 



m 



s(l) 



(27r)3 



q(q-k) 



Iq-kr |q-k 

(B13) 
(B14) 



Further progress is made by noting that it is not 
necessary to give a detailed account of the bilinear 
mode-coupling term in m'^^°(f) if one is interested in 
the anomalous dynamics near the diffusion-localization 
threshold. Indeed, on the one hand, it is clear that this 
term can be singular when rif is very small, becaiise the 
diffusion-localization and liquid-glass transition lines are 
then very close, but it is also very small, because of the 

. (2) 

explicit factor Uf in Wj, . In particular, it vanishes iden- 
tically when nf = 0. It can thus be neglected for small 
nf. On the other hand, for larger nf, but not near point 
T which would require a refined treatment, the wide sep- 
aration between the two transition lines guarantees that 
the coll^tive dynamics will be fast and regular at the 
PP' cj|^'^|iAJlocalization transition. It results that the sin- 
''^ular laKed-particle dynamics is cut off in the bilinear 
term, wmch generates a fast and smoothly varying con- 
tribution to m^^^{t). Thus it is in general enough to 
set 



k'^dk £^ss(b) 



\it) (C4) 



The c'g^^'^^'s and /iq'^''^ are the Fourier transforms of 
the connected single-particle-fluid direct correlation func- 
tions and of the single-particle-single-particle blocked to- 
tal correlation function, respectively [32] . 

The equations reported in the main body of the paper 
are immediately obtained for n — 1. 



APPENDIX C: HYDRODYNAMIC 
APPROXIMATION FOR THE 
DIFFUSION-LOCALIZATION TRANSITION 

In this appendix, the results of a generalized hydrody- 
namic approximation, in which the wave vector depen- 
dence of the MCT equations is simplified, are reported for 
the dynamics in the vicinity of the diffusion-localization 
transition. They represent an extension of Leutheusser's 
theory [26] to finite fluid densities. 

The theory is best obtained by considering Laplace 
transformed dynamical quantities, defined as (p{z) = 
i dte"*(j){t) for Im ^; > 0. Equation (11) then be- 
comes 



with an appropriate renormalization of 7. Then, combin- 
ing Eqs. (CI) with Dq{z) ~ D{z), (C3), and the Laplace 
transform of Eq. (C4), a transcendental equation for D{z) 
is obtained, which reads 



z + ij^ b{z) = ~1 + 



+°° k^dk k 



ss(b) 



67r2 k^ + z/D{z)' 
(C5) 

Prom the solution of this equation, which is the core re- 
sult of the proposed generalized hydrodynamic approx- 
imation, all relevant dynamical quantities can be com- 
puted. 

The diffusive and localized phases are characterized by 
different small z behaviors of D{z). In the diffusive phase, 



lim D{z) 

z—>iO 



iD, 



(C6) 



where D is the finite long-time diffusion coefficient. Prom 
Eq. (C5), one finds 



D = 



1 



i/i««(b)(r. = 0) 



-1 



q^Dgiz) 



(CI) 



In the localized phase, 



D{z) ~ zrf, 



(C7) 



(C8) 
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where n is the finite localization length, which obeys 



p + ao 

Jo 



k*dk h] 



s{h) 



67r2 A;2+rf2 



(C9) 



So, as anticipated, the present theory indeed describes a 
bona fide difFusion-locahzation transition scenario, with 
the above condition providing an implicit equation of the 
transition line. 



Using a large ri expansion of Eq. (C9), one easily shows 
that the divergence of r\ and the vanishing of D both 
occur when 



/i=^(b)(r = 0) = 3. 
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